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Abstract 

We present analytic formulae for the QCD renormalization group factors relating the Wilson 
coefficients Ci{iit) and Ci{fj,), with /i^ = 0{mt) and ji < fit, of the AF = 2 dimension six four- 
quark operators Qi in the Standard Model and in all of its extensions. Analogous analytic 
formulae for the QCD factors relating the matrix elements {Qi{2 GeV)) and {Qiifix)) with 
/Ik < 2 GeV are also presented. The formulae are given in the NDR scheme. The strongest 
renormalization-group effects are found for the operators with the Dirac structures (1 — 75) ® 
(1+75) and (1— 75)®(1— 75). We calculate the matrix elements {K^\Qi\K^) in the NDR scheme 
using the lattice results in the LRI scheme. We give expressions for the mass differences AMk 
and AMb and the CP-violating parameter ex in terms of the non-perturbative parameters 
Bi and the Wilson coefficients Ci{fit)- The latter summarize the dependence on new physics 
contributions. 



1 Introduction 



Renormalization group short-distance QCD effects play an important role in — and 
B^g — B^ij^g mixing within the Standard Model (SM) and its extensions |jl], |2|. They can be 
calculated by solving renormalization group equations that govern the scale dependence of the 
Wilson coefficients Cj(yu) of the relevant AF = 2 operators Qi. The resulting effective weak 
Hamiltonian reads 

H^i=' = ^,M'w E ^CKMa(/i)Q. . (1.1) 

Here is the Fermi constant and V^j^^^ the Cabibbo-Kobayashi-Mask;awa (CKM) factor equal 
to {y^^VtdY the case of — fi^ mixing in the SM. Beyond the SM other factors not pro- 
portional to CKM elements are generally present. Using this Hamiltonian one can calculate 
AF = 2 amplitudes, in particular the mass differences AMk and AM^^^ in the — and 
B^^ — B^^g systems and the CP- violating parameter sk- 
Within the SM there is only one single operator 

gVLL ^ (s^^^p^ci") (s^^Pirf^) (1.2) 

relevant for — mixing, with analogous operators for B^^ — ^ mixing obtained from 
(|1.2|) through the appropriate change of flavours. Beyond the SM the full set of dimension 
six operators contributing to — mixing consists of 8 operators that can be split into 5 
separate sectors according to the chirality of the quark fields they contain. These operators are 
listed in (|2.1|). Corresponding operators contributing to B^g — 5^,. mixing exist. 
The general expression for Ci(/i) is given by 

c{|l) = iJ{^l,^lt)C{^lt) (i.s) 

— * ^ 

where C is a column vector built out of the Cj's and U{iJ,,Ht) is the renormalization group 
matrix. C{fit), with fit = 0{mt), are the initial conditions which depend on the short distance 
physics at high energy scales. In particular they depend on the top quark mass and the 
couplings and masses of new particles in extensions of the SM. We will later briefly discuss the 



case of scales much higher than rrit. Otherwise Ht denotes a high energy scale in the range, say, 
Mw < < 2mt. 

While the initial conditions Ci{fit) at the NLO level are known only in the SM 0, ^ and in 
some of its extensions 0], all the ingredients are available to compute the NLO evolution matrix 
U {fi, fit) for all possible extensions of the SM. Indeed, the two-loop anomalous dimension matrix 
for all AF = 2 four-quark dimension six operators has been calculated in the regularization 
independent renormalization scheme (RI) in |^ and in the NDR scheme in P]. Together with 
the known one-loop anomalous dimension matrix 0, § and the known (3 function, the evolution 
matrix can be straightforwardly computed by means of the methods reviewed in 0. 

The LO analytic expressions for f/(/i, /i^) can be found in ^j. For phenomenological appli- 
cations it is useful to derive analogous expressions including NLO corrections. The first step in 
this direction has been made in where U {fi, fig) with fis > rrit has been written as 

U {fi, fis) = U {fi, fit)U {fit, fis) (1.4) 

with U {fit, fis) given analytically in the Landau RI scheme (LRI) but (j{fi,fit) evaluated nu- 
merically for fi = 2 GeV, fit = frtt and particular values of iric, mt and as- The corresponding 
results for i?^^ — _B°^ mixing have not been presented in 0. 

It should be emphasized that NLO corrections are necessary for a satisfactory matching of 
the Wilson coefficients to the matrix elements obtained from lattice calculations. Moreover 
as demonstrated in p|, ^ the inclusion of QCD corrections at the LO and the NLO level is 
mandatory in order to place reliable constraints on the parameters in the extensions of the SM, 
in particular on the squark mass matrices in supersymmetric theories. 

The purpose of our paper is to present NLO analytic formulae for the matrix U {fi, fit) relevant 
for B^g — 5^] g mixing {fi = 0{mi,)), and — mixing {fi = 0{1 — 2 GeV)). These formulae 
when combined with the initial conditions C{fit) and the hadronic matrix elements {Q{fi)) will 
allow to calculate in the future the AF = 2 amplitudes for any extension of the SM. 

The formulae given below for U {fi, fit) apply to the situation in which the initial conditions 
for the Wilson coefficients are known at fit = 0{mt) and the evolution down to scales fi < fit 
is performed in an effective theory with the top quark and the heavy new particles integrated 
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out. Whether the top quark and the new particles have been integrated out at a single scale 
fit or at different scales, say /if, /is^, /isa with fit < jis^ < jis-^i is immaterial here. What matters 
are the values of the Wilson coefficients at /it and not how they have been evaluated from the 
contributions at scales higher than /i^. On the other hand in the process of the evaluation 
of Ci^jit) large logarithms log fis^/ fit, log/isa/z^si may appear. These logarithms have to be 
resummed which results in new evolution functions f/(/ii,/isJ, ?7(/is^, /Xsj), etc. As discussed in 
[H, 1^ the structure of these matrices is model dependent and consequently beyond the scope 
of the present paper. We will, however, provide an analytic formula for the evolution U{^t, /^s) 
with /it -C fis in an effective f = Q theory in which only SM degrees of freedom are present and 
all new particles have been integrated out. 

Now, the lattice results for the matrix elements {Q{fi)) are usually given at fi = 2 GeV. In 
what follows we will denote this scale by fii. On the other hand large- approaches, the chiral 
quark model and any non-perturbative method in which the low-energy degrees of freedom are 
mesons provide these matrix elements at scales fix < IGeV. In our opinion it would be useful 
to have the matrix elements obtained by means of different methods at a common "standard" 
scale, which we will choose to be fi^ in the following. This is achieved using the formula 

{Q{fiL)) = U''if^K,l^L){QM) (1.5) 

where {Q{^k)) are the matrix elements calculated for scales < /^l and (J{^k-,^'l) is 
the renormalization group evolution matrix. In our paper we provide analytic formulae for 

At this point we would like to stress that our paper is addressed first of all to the practitioners 
of weak decays who do not want to get involved with the details of NLO calculations but 
rather would like to use the final QCD factors in phenomenological applications. On the other 
hand it should also be useful to experts. Indeed, having explicit analytic formulae, rather than 
numerical values, not only gives the freedom to change input parameters but also makes possible 



the checking of a given calculation. In particular when multiplying the U matrices like in (1.4) 
one easily generates higher-order terms in as which really do not belong to NLO corrections. 
While these corrections should be removed from NLO expressions, this is not always done in 



the literature. Consequently, already at this stage unnecessary discrepancies of the order of 
5% between calculations performed by different groups may arise. These higher-order terms in 
a<j are consistently removed in the present paper. We are aware of the fact that some of the 
formulae presented below are rather long. Nevertheless we believe that they should turn out 
to be useful in future phenomenological applications. 

The paper is organized as follows. In Section 2 we give the list of the AF = 2 operators 
in question and establish our notation. In Section 3 we give analytic formulae for the QCD 
factors [f]ij{f^)]a that represent the evolution matrix U{^,fit) in ( |1.3|) in five different sectors, 
a = (VLL, LR, SLL, VRR, SRR), in the leading order (LO) and the next-to-leading (NLO) 
approximation in the NDR scheme. In Section 4 we give the analogous formulae for the QCD 
factors [Pij{fJ'K)]a which represent the evolution matrix U^hk^IJ'l) in (|1.5| ). In Section 5 we 
provide numerical results for [?7ij(At)]^ and [pjj(yu)]^ in the NDR scheme. In section 6 we dis- 
cuss the transformation rules for obtaining the corresponding results in other renormalization 
schemes and we present the relation between the QCD factors calculated here and the QCD 
factors TjB and ri2 used in phenomenological applications. In Section 7 we calculate the matrix 
elements {l^\Qi\K^) in the NDR scheme using the lattice results in the LRI scheme j^, p!0|] . 
We give general expressions for the mass differences AM^ and AMb and the CP-violating pa- 
rameter ex in terms of the non-perturbative parameters and the Wilson coefficients Ci{fit)- 
We conclude in Section 8. For completeness we list in appendix ^ the one-loop and two- loop 
anomalous dimension matrices that we have used in our paper. Appendix ^ gives the general 
formulae for the U matrices which have been used to obtain the analytic formulae of sections ^ 
and ^ Finally in Appendix C we give analytic formulae for the evolution matrix U{pt,fis). 

2 Basic Formulae 

For definiteness, we will give explicit expressions for the operators responsible for the K'^ — 
mixing. The operators belonging to the VLL, LR and SLL sectors read 
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LR 



Q 



LR 



SLL 



Q 



SLL 



{rPLd''){s^PRd^), 

{rPLd''){s^PLd'^), 
ira,,PLd''){s''a>'''PLd''), 



(2.1) 



where a, /3 are colour indices, cr^,y = |[7/^,7v] and P^^/j = | (1^=75). The operators belonging to 
the two remaining sectors (VRR and SRR) are obtained from Q^^^ and Qf^^ by interchang- 
ing Pl and Pr. Since QCD preserves chirality, there is no mixing between different sectors. 
Moreover, the anomalous dimension matrices and the evolution matrices in the VRR and SRR 
sectors are the same as in the VLL and SLL sectors, respectively. Therefore, in the following, 
we shall consider only the VLL, LR and SLL sectors. However, one should remember that the 
initial conditions Ci{^t) are generally changed when P^ and Pr are interchanged. The oper- 
ators in the case of P° — p|j mixing are obtained from (|2.1|) through the replacement s ^ h. 
Performing the subsequent replacement d s gives the operators contributing to P^ — P^ 



mixing. The one-loop and two-loop anomalous dimension matrices of the operators (|2.1|) are 
given in appendix 

Restricting the discussion to the VLL, LR and SLL sectors f/(/ii,/i2) takes the following 
form 

\ 



( [^(/^l,/^2)]vLL 

[^(/^i,/^2)]lr 

V [^(/^l./^2)]sLL 



(2.2) 



where [r)(/ii, /i2)]LR and [57(Ati, /U2)]sll 2x2 matrices and /ii < /i2- In what follows we will 
use a short-hand notation, denoting the QCD factors representing f/(/i, /x^) and U{fiK,fJ'L) by 



[p{f^K,f^L)]a 



+ 



c^f\p 



Arc 



r)«(/i) 



(2.3) 
(2.4) 



respectively. That is, we will suppress the high-energy scale fit in the argument of the r^-factors. 
Similarly, we will suppress the "lattice scale" /xl in the argument of the p-factors. Using this 



notation we have for instance 



VLL/ 



[^ll(^b)]LR iVuif^b)] LR 
[^2l(^b)]LR [^22(/ife)] LR 



(QrM) = [p(/ix)]vLL(Qr(/^i^)) 



LR 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



V [Pl2(Pi^)]LR [P22(/ii^)] LR , 

and analogous formulae for the SLL sector. Note that in accordance with (|1.5| ), the transpose 
of [p(pi<')]LR enters the transformation ( |2.8| ). 

In Section 3 we will give analytic formulae for the LO factors r]^^\fi) 



vll' 



lr' 



SLL 



and the NLO factors 



r^«(/i) 



VLL 



LR 



SLL 



The corresponding 



expressions for the p-factors are given in Section 4. a^/\fi) is the QCD coupling constant in an 
effective theory with / flavours: / = 5 for /it < /i < /i^, / = 4 for /ic < p < Pfe, and / = 3 for 
fi < ^c, where fit = 0{mt), fib = 0{mb) and fic = O^rric)- We impose the continuity relations 



(2.9) 



The general expression for a^/^ reads 



«F(P) 
47r 



1 



/3oln U2/A. 



MS 



with 



11-3/, 



38 



In^ U2/A. 



MS 



(2.10) 



A = 102 - -/, 



is the QCD scale parameter in a theory with / quark flavours |11 



(2.11) 



The exist- 



A. 



(5^ 
MS 



(226 )MeV. 



ing analyses of high energy processes give a-^(M^) = 0.118 ± 0.003 [0 or equivalently 



The evolution matrix, U {fx, fit), is given as follows: 
U{fi,fit)=T,e.v l^^/9j^ 



(2.12) 
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with g denoting the QCD effective couphng constant and Tg an ordering operation defined in 
I3{g) governs the evolution of g and 7 is the anomalous dimension matrix of the operators 
involved. 

We also have 

IJ{HK,IJ't) = U{fXK, fJ'cWiHc, flb)U{flb, fJ't) (2.13) 

with the three factors on the r.h.s. evaluated in / = 3, / = 4 and / = 5 effective theories, 
respectively. Now, 

U {fiL, fJ't) = U ifJ'L, Hb)U{fib, fJ't), U ifJ'K, fJ't) = U (/iK, fiL)U{fJ'L, /it). (2.14) 

This means that knowing rjij^fiL) and Pij^fix) allows to calculate rjij^fix)- 

Keeping the first two terms in the expansions of ^{g) and j3{g) in powers of g 

inserting these expansions into ( |2.12| ) and ( |2.13| ) and expanding in as one can calculate the 1]- 



and p-factors defined in (|2.3|) and ( p. 41) , respectively. To this end, one has to remove terms 
(9(a;^) and higher-order terms. We discuss this point in appendix where the expansions of 
U (/i, fit) for = pb, fJ' = fi'L and p = are given. 

3 77-Factors in the NDR Scheme 
3.1 ?7-Factors for B^^ — mixing 

VLL-Sector 



= (3.1) 



6/23 



VLL 



^^'^^^''^JvLL = 1-6273(1 -r^5)r]r'. (3.2) 



LR-Sector 



^n{^^b)\^^ = 41 '\ (3.3) 
'nf2{^b)\^^ = 0, (3.4) 



3/23 





?7oi^ (Uh) 


LR 


= 2. 3/23 _ -24/23. 

^ \ 10 10 J ' 


(3.5) 






LR 


-24/23 
= ^75 ) 


(3.6) 






LR 


= 0.9250775-2^/2^ + r?5^/23 ('_2.0994 + 1.1744%) , 


(3.7) 

V / 




^12 Vb) 


T R 

Ijri, 


= 1.3875 (r^s''/'' - r/s'^^/^s) , 


(3.8) 






T R 


= (-11.7329 + 0.7829 r^s) r/g^/^s ^ ^^.-24/23 (_5,3048 + 16.2548 r^s) , 


(3.9) 






LR 


= (7.9572 - 8.8822 r/s) r^g-^Vss ^ 0.92507752^/23. 


(3.10) 


SLL-Sector 












SLL 


= 1.0153r/5-o-63i5 _ o.0153r75°•^l8^ 


(3.11) 






SLL 


= 1.9325(r75-°-6='i5- 770-^184)^ 


(3.12) 




^21 Vb) 


SLL 


= 0.0081(r/°-^i«^-r75-°-63i5), 


(3.13) 






SLL 


= 1.0153r/°-^^^^ - 0.0153^75-°•6=^^^ 


(3.14) 




Vll{^^b) 


SLL 


= (4.8177 - 5.2272 r/5) 775-°-^=^^^ + (0.3371 + 0.0724 r]^) r/5°•^^^^ 


(3.15) 




Vuif^b) 


SLL 


= (9.1696 - 38.8778 rj^) rj^^^-^^^^ + (42.5021 - 12.7939 r/5) ri5^-™\ 


(3.16) 




V2l{f^b) 


SLL 


= (0.0531 + 0.0415 r/5) r75-° *^3i5 _ (o.0566 + 0.0380 r/5) ^75°•^^^^ 


(3.17) 




^22 Vb) 


SLL 


= (0.1011 + 0.3083 ri5) rj^-^-^^^^ + 775°-^^®^ (-7.1314 + 6.7219 r/5) , 


(3.18) 



where 



V5 



(3.19) 



3.2 77- Factors for K — K mixing with fi = fii 



These factors are relevant for — Tf mixing but can also be used in D° — mixing. They 
can be expressed in terms of rj^ defined in (|3.19|) and 



r]4 = 



VLL-Sector 



VLL 
VLL 



6/25 6/23 

ril^"^^ 4^^^1.7917 - O.WAArii - 1. 6273774^75) • 



(3.20) 

(3.21) 
(3.22) 



LR-Sector 



'111 Kr'J^J 


LR 


3/25 3/23 
'/4 '15 -1 


(3.23) 


ni2} ( ij,T^ 
'/12 yf^L; 


LR 


= 


(3 24) 


V2I [f^L) 


LR 


2 / 3/25 3/23 -24/25 -24/23n 

1 VI ' VI / VI ' , VI ' 1 

— -^yli V5 ^ Vi V5 ): 


yO.Z,o ) 


'122 KH'L) 


LR 


-24/25 -24/23 
'/4 'lb I 


f3 261 


'111 Kh-l) 


LR 


- n 9279 riA -^^/^^ 73. -2^/23 _ A QQoq r?. ^^/^^ 77. "^^/^^ 








\n 3/25 3/23 / 9 ri9/( l ri nyc^'? n _L 1 17AA n n \ 




'112 KH'L) 


LR 


- _1 "^qi fi -24/25 „ -24/23 , n nn^o „ 28/25 „ -24/23 








+1.3875 r/4^/^^ ?75 , 




V21 (i^l) 


LR 


= -0.0019 7^4''/'' r?5-''/'' + 5.0000 774^/'' r^s'/'' 


(3.29) 






+7?4^/^^ 775^/2^ (-16.6828 - 0.0502 774 + 0.7829 774 775) 








+774"^^/^^ 775"^^/^^ (-4.4701 - 0.8327774 + 16.2548 7/4 775) , 




V22 (i^l) 


LR 


- 0.0029 774^^/2^ 775-2^/23 + 0.9250 7742^/2^ 775^^/2^ 


(3.30) 






+774"^^/^^ 775-^^/^^ (6.7052 + 1.2491 774 - 8.8822 774 775) . 





SLL-Sector 





SLL 


= 1.01537/4-°-^«i'^775-°-63i5 - 0.015377°•6«l°77°■^l«^ 


(3.31) 


Vl2{f^L) 


SLL 


= 1.9325(7/4-o-58i%5-°-«=^i^ -7;4°-66io^o.n84)^ 


(3.32) 


V21W) 


SLL 


= 0.0081(77°-66^°770-^^«^ - 774-°-5«^V°-''''), 


(3.33) 




SLL 


= imf>H'''V,''^ - 0.0 1 537;4-°-5«^V°'''', 


(3.34) 




SLL 


= 0.0020 774^-^^^° 775-°-*^=^^^ - 0.0334 774°-^i90 ^^0.7184 


(3.35) 






+774-°-^^^° 775-°-^^^^ (4.2458 + 0.5700 774 - 5.2272 774 775) 








^^^0.6610^^0.7184 (0.3640 + 0.00647/4 + 0.07247/47/5) , 




^12 Vl) 


SLL 


= 0.0038774^-^'^° 775-°-''^' - 4.2075 774°-'i9° 775°-'''' 


(3.36) 






^^^-0.5810^^-0.6315 (g_0810 + 1.08487/4 - 38.87787/4 7/5) 








+7/4°-^''^° 7/5°'^^^^ (45.9008 + 0.80877/4 - 12.7939 7/4 7/5) , 
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SLL 



-0.0011 774'-'''° r75"°-'''' + 0.0003 774°'^'° r/5°-'''' 
+r/4°-66^° r?5°-^^^^ (-0.0534 - 0.0034 774 - 0.0380 r/4 r^g) 
^^^-0.5810 ^^-0.6315 (0.0587 - 0.0045 r]^ + 0.0415 r/4 r/5) 



SLL 



-0.0020 7/4^-^^^° r/5-°-63^5 



+ 0.0334 V-'i9%5°-'''' 



^^^-0.5810 ^^-0.6315 (0.1117 - 0.0086 ?74 + 0.3083 r]^ r]^) 
+r/4°-^'^^° 775°-^^^^ (-6.7398 - 0.4249 7/4 + 6.7219 r]^ 775) 



(3.37) 



(3.38) 



3.3 ?7-Factors for — K mixing with ^ik = C(lGeV) 



The formulae for the QCD factors fj^fix, fit) = V^I^k) relating the coefficients Ci^fix) and Ci{fit) 
are rather long and will not be presented. These factors can be obtained using the relation 



V{1^k) = p{flK)fl{flL) 



(3.39) 



with fj{pi) given in section and pif^x) given below. When calculating ( p. 391) , terms of 0(«g) 
should be removed. 



4 p-Factors in the NDR Scheme 



These factors allow to calculate {Qi^Ph)) from {Qi^Pk)) with px < Pc- They can be expressed 
in terms of 



= 



VLL-Sector 



and 



af\pc 



(4.1) 



P^'\pk) 



VLL 
VLL 



2/9-6/25 
^3 Va , 

1.8951r^f r^4^/^^ - ^ .lQ?>?>r,T f,T ' l.lQnvf'vf"'- 



(4.2) 
(4.3) 



LR-Sector 



P12V 



K 



LR 



LR 



1/9-3/25 



(4.4) 
(4.5) 
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LR 
LR 
LR 

LR 

LR 



LR 



SLL-Sector 



2/ 1/9-3/25 
-8/9 ~-24/25 



8/9 — 24/25 



0.9306 773-'/' 774-''/'' - 0.0027 r?3''^/'' 774~''^'' 
+m^'^fiT^ (-1.9784 - 0.0457r/3 + 1.0962 773 , 
-1.3958 7^3-'/' f,f"^' + 0.0040 r^,'"' f,f"^' 

+1.3918 7^3'°/' ^74'^'', 
-3.8570 r/3-'/' 7^4-''/'' 
+r?3^/' r^^^^/^^ (-0.6113 - 0.0018 773 + 20.4220 774) 
+??3^/' r/f (-16.6523 - 0.7407 log(773) - 0.0305 773 + 0.7308 773 774) 



(4.6) 
(4.7) 
(4.8) 

(4.9) 

(4.10) 



5.7855 773-'/' 774-'^/'^ + 0.9279 773^°/' 774''/'' 

fjf^^^^ (0.9170 + 0.0027773 - 7.6331 774) 



(4.11) 



Pii{Pk) 


SLL 


= 1.0153773-°-53^'774-°-^«^° - 0.0153773°-«^^°774°-^^^°, 


(4.12) 


PuiPK) 


SLL 


= 1.9325(773-°-5379~-0.5810 _ ^0.6120-0.6610)^ 


(4.13) 


P2i{Pk) 


SLL 


= 0.0081 (773°-6120 -0.6610 _ ^-0.5379 --0.5810) ^ 


(4.14) 


P22{Pk) 


SLL 


= 1.0153773°-«i2'^77°-6«i° - 0.0153773-°-5379~-o.58io^ 


(4.15) 


PiKpk) 


SLL 


= 0.0019 773^-^1^° 774-°-5«^° - 0.0663 773°-'^'^ vl'''' 


(4.16) 






^^^-0.5379 ^j0.58io (3 g4g7 + 0.3952 773 - 4.6906773 774) 








+^30-6120 -0.6610 (0.4264 + 0.0039 773 + 0.0808 773 774) , 




PuiPK) 


SLL 


= 0.0036 773^-^^^° 774-°-^«i° - 8.3647 773°-^^^^ ^"'^^^^ 


(4.17) 






+,^3-0-5379 j^jo.5810 (7 3253 + 0.7521 773 - 42.0005 773 774) 








+773°-^^^° 775-^^^° (53.7722 + 0.4933773 - 11.9813 773 774) , 




P2i{Pk) 


SLL 


= -0.0010 773^-^^'° 774-°-5si° + 0.0005 773°-^6'^ 71^'''' 


(4.18) 






^^^0.6120 -0.6610 (_o.0519 - 0.0021 773 - 0.0425 773 774) 








+^3-0-5379 ^^-0.5810 (0.0628 - 0.0031 773 + 0.0372 773 774) , 
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P^k^K)^, = -0mi9v3'-''''v^'-'''' + 0M63v3'''''vl-'''' (4.19) 
^^^-0.5379 ^j0.58io (^0.1196 - 0.0060 ris + 0.3331 r^s r/4) 

+ri3°-^^^^ (-6.5470 - 0.2592 r/3 + 6.2950 r/3 Vi) ■ 



Regarding the appearance of log (r^s) in eq. ([4.10|) in the LR sector, we direct the reader's 



attention to eq. (|B.4| ) and the fact that in the LR sector for / = 3 the form of the evolution 
operator given there breaks down, as the matrix J has a singularity at / = 3. However, taking 
the limit / — > 3 of the complete expression ( |B.4| ), a finite result exhibiting the aforementioned 
term 0{as) log (r^s) is obtained [jl3|. In accordance with the convention of ( p.6|) , ( p.7| ) we treat 
it as an NLO contribution to the evolution operator. 

5 Numerical Results 

In tables |1|-^ we give the numerical values for the rjij and pij factors in the NDR scheme. To 
this end we have used a^^\Mz) = 0.118 ± 0.003 with the corresponding values of a^/^ and 
A^g^ for f = 4: and / = 3 theories, see table 2 in the first paper in 0. Moreover we have set 
fit = 'fntij^t) = 166 GeV, /it = 4.4 GeV, ji^ = 2.0 GeV, /ic = 1.3 GeV and fix = 1-0 GeV. In 
order to illustrate the effect of the NLO corrections we show also the results in the LO. In doing 
this we have, however, used the two-loop expression for in both the LO and the NLO parts. 
In figs. |l|andHwe show the factors [f]ij{p)]i,B. and [flijip)]sLL versus /i setting a^^\Mz) = 0.118. 

Let us recall that in the absence of renormalization group effects [t]{p)]yll ~ [p(/^)]vll ~ ^ 
and [flifJ')]^ and [p(yu)]^ are unit matrices. Renormalization group effects generate non-diagonal 
elements in these matrices and renormalize ['7(yu)]vLL and the diagonal elements in [i7(/i)]„ and 
[p{p)]a away from unity. 

Inspecting tables and figs. we observe the following pattern: 

• Large renormalization group effects are found for the diagonal entries [?722(/^)]lr and 
[^ii(/^)]sLL "which for af'^Mz) = 0.118 and yU = hk = 1 GeV are enhanced by factors of 
5.6 and 2.9, respectively. On the other hand ['?22(Ai)]sLL strongly suppressed down to 
0.26 at fiK = I GeV. 
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• Similarly the renormalization group effects in the non-diagonal entries [?72i(/^)]lr and 
[^i2(Ai)]sLL large. For fix = 1-0 GeV and a'^^\Mz) = 0.118 they become as large as 
—3.2 and 4.8, respectively. This implies that C2^(/x) is strongly affected by the presence 
of the operator Q\^. Similarly Cf^^(/x) is strongly affected by the presence of 

• These enhancements and suppressions are more pronounced after the inclusion of NLO 
corrections. The largest NLO corrections, in the ball park of 25%, are found for the 
elements [?72i(/^k)]lr, [?722(//k)]lr and [?722(/^k)]sll- 

• [^7ii(a*)]lr and ?7vll(a*) are both suppressed but this suppression is at most by 10% and 
30%, respectively. 

• [^i2(Ai)]LR is small and [?72i(a*)]sll negligible in the full range of n considered. This implies 
that Ci^{ijl) and C|^^(/x) are essentially unaffected by the presence of the operators 
and Qf^^, respectively. 

• Similar patterns are observed for the factors. 

On the basis of this pattern we conclude that the renormalization group effects strongly enhance 
the Wilson coefficients 6*2^ and Cf^^ and strongly suppress C|^^ with respect to their values 
at fjLf The corresponding effects in CY^^ and are substantially smaller. 
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a(5)(M2) = 0.115 


af\Mz) =0.118 


af\Mz) = 0.121 




LO 


NLO 


LO 


NLO 


LO 


NLO 


[^(/^fe)]vTT 

L'\' ^ / J \ 1 ; 1 ; 


0.835 


0.847 


0.829 


0.842 


0.823 


0.836 


['7ll(/ife)]TR 
L "/J 1 J 1 \ 


0.914 


0.923 


0.911 


0.921 


0.907 


0.919 


[^12(/^fe)]LR 





-0.037 





-0.041 





-0.045 


[^2l(/^b)]LR 


-0.760 


-0.835 


-0.801 


-0.885 


-0.845 


-0.939 


[?722(/^b)]LR 


2.054 


2.181 


2.112 


2.254 


2.176 


2.334 


[^ll(/^b)]sLL 


1.560 


1.621 


1.587 


1.654 


1.616 


1.690 


[^12(/^6)]SLL 


1.809 


1.910 


1.884 


1.993 


1.962 


2.082 


[^2l(/^6)]sLL 


-0.008 


-0.006 


-0.008 


-0.007 


-0.008 


-0.007 


[r/22(/^6)]sLL 


0.595 


0.563 


0.583 


0.549 


0.570 


0.535 



Table 1: Numerical values for the 77-factors for i?^^ — i?^^ mixing. 







) = 0.115 


af^{Mz 


) = 0.118 


cSp{Mz 


) = 0.121 




LO 


NLO 


LO 


NLO 


LO 


NLO 


[^(/^l)]vll 


0.778 


0.796 


0.768 


0.788 


0.757 


0.780 


[?7ii(/^l)]lr 


0.882 


0.906 


0.876 


0.906 


0.870 


0.906 


[7712(/^l)]lR 





-0.076 





-0.087 





-0.101 


[r72l(//L)]LR 


-1.236 


-1.401 


-1.336 


-1.530 


-1.449 


-1.677 




2.735 


3.009 


2.879 


3.200 


3.043 


3.420 


[^ll(^i)]sLL 


1.859 


1.976 


1.918 


2.052 


1.984 


2.138 


[m2(/^L)]sLL 


2.586 


2.775 


2.732 


2.946 


2.892 


3.137 


[r?2l(/^L)]sLL 


-0.011 


-0.009 


-0.011 


-0.009 


-0.012 


-0.009 


[?722(/^l)]slL 


0.480 


0.438 


0.461 


0.417 


0.442 


0.394 



Table 2: Numerical values for the r;-factors for — K mixing with jji = = 2 GeV. 
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) = 0.115 


a(5)(Mz) = 0.118 




) = 0.121 




LO 


NLO 


LO 


NLO 


LO 


NLO 


L / Vr^JA- / J VLL 


0.701 


0.735 


0.681 


0.720 


0.658 


0.705 


[^711 (Uk)]j-d 
I ILL \f^j\ yjLH 


0.837 


0.921 


0.825 


0.941 


0.811 


0.978 







-0.194 





-0.254 





-0.347 


[?72i(/^k)]lr 


-2.199 


-2.657 


-2.545 


-3.159 


-3.006 


-3.861 


[?722(/^x)]lR 


4.136 


4.875 


4.643 


5.630 


5.320 


6.688 


[r?ll(/^K)]sLL 


2.392 


2.663 


2.566 


2.912 


2.787 


3.243 


[^12(/ii^)]gLL 


3.836 


4.273 


4.223 


4.782 


4.702 


5.442 


[^2l(/^i^)]gLL 


-0.016 


-0.011 


-0.018 


-0.012 


-0.020 


-0.012 


[r/22(/^i^)]sLL 


0.346 


0.291 


0.314 


0.255 


0.279 


0.217 



Table 3: Numerical values for the 77-factors for — K mixing with ijLk — 1 GeV. 





afKMz 


) = 0.115 




) = 0.118 




) = 0.121 




LO 


NLO 


LO 


NLO 


LO 


NLO 


[p(/^i^)]vLL 


0.902 


0.923 


0.887 


0.914 


0.870 


0.905 




0.950 


0.955 


0.942 


0.951 


0.933 


0.947 


[pi2(//if)]LR 





-0.045 





-0.060 





-0.083 


[P2i(/^x)]lr 


-0.375 


-0.448 


-0.447 


-0.548 


-0.544 


-0.688 


[P22(^x)]lR 


1.512 


1.620 


1.613 


1.762 


1.748 


1.963 


[pil(/^x)]sLL 


1.293 


1.357 


1.345 


1.431 


1.413 


1.533 


[pi2(/ix)]sLL 


1.029 


1.176 


1.190 


1.381 


1.394 


1.650 


[p2l(/^i^)]sLL 


-0.004 


-0.003 


-0.005 


-0.003 


-0.006 


-0.003 


[P22i^iK)]sLL 


0.744 


0.688 


0.710 


0.641 


0.670 


0.584 



Table 4: Numerical values for the p-factors with hk = 1 GeV and /x^ = 2 GeV. 
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Figure 1: The [%]lr factors as functions of // for af\Mz) =0.118. 
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Figure 2: The [%]gLL factors as functions of /i for af \Mz) = 0.118 
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6 General Remarks 



6.1 Renormalization Scheme Dependence 

The evolution matrix U is renormalization scheme dependent. It is instructive to recall how 
this scheme dependence is canceled in physical amplitudes by considering a single operator Q. 
Then the AF = 2 amplitude reads 



A{AF = 2) = {H,s) = ^^MI,VckmC{^){Q{^)) . 



(6.1) 



The Wilson coefficient is given by 



c{^^) = u{^^,^^t)C{^lt) 



(6.2) 



where 



1 + ^ J 

4:71 



(6.3) 



with 



7 



(0) 



Wo 



P 7(1) 



(6.4) 



and 



(6.5) 



Cq and Ci depend generally on m^, M]y and the masses of new particles in extensions of the 
SM. 

Now, the renormalization scheme dependence of Ci is canceled by the one of J in the last 
square bracket in (|6.3| ). The scheme dependence of J in the first square bracket in (|6.3|) is 
canceled by the scheme dependence of {Q{fi)). The power P and the coefficient Co are scheme 
independent. 



6.2 Transformation to Different Renormalization Schemes 



Once the Wilson coefficients Cj(yu) have been calculated in the NDR scheme, they can be 
transformed to a different renormalization scheme A by means of 



1 



An 



NDR^A 



(6.6) 



Likewise the matrix elements {Qi{fi)) can be transformed from scheme A to the NDR scheme: 



NDR 



+ 7'"' AfA^NDR (<5(yU))A , AfA^NDR = -AfNDR-.A • 

An / 



(6.7) 



The transformation matrices ArNOR^Ri from the NDR scheme to the RI schemes of can be 
found in section 5 of UBI. 



6.3 r]B and 772 Factors in the SM 

At this point we would like to warn the reader that the QCD factors rjB = 0.55 P, |^ and ri2 
used in the analysis of B'^ g — mixing and of the top quark contribution to Ek in the SM 
should not be identified with the factors [?7(yUfe)]vLL ['7(/^i^)]vLL presented in this paper. 

The factors rjB and 772 are discussed in detail in [||, ||, |^. See in particular the expressions 
(12.10) and (13.3) in [||] for 772 and rjB, respectively. Using these expressions it is straightforward 
to find the relation between rjB and [v{fJ'b)]yi^i^- It reads 



VLL '-'SM 



-6/23 



rjBASoixt) 



(6.8) 



where C^u^ifit) includes NLO corrections calculated in [0, 3, J5 = 1.627 in the NDR scheme 
and ASo{xt) with xt = m^(/ii)/M^ is the LO expression for Cg^^(/it) that is obtained from box 
diagrams with top quark exchanges without QCD corrections. rjB in contrast to [^7(Ai6)]vLL 
renormalization scheme independent and does not depend on /i^. The latter dependence has 
been factored out as seen on the r.h.s of (|6.8|). Notice that the QCD corrections to C^u^ifit) 
have been absorbed into rjB. An analogous relation between 772 and [^7(Aii^)]vLL t)e found. 
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6.4 Going Beyond the SM 

In the SM 



where Bb is the renormahzation group invariant parameter defined by 



5b = S^LL(^,) «(5)(^, 



'6/23 



1 + 



An 



-J. 



(6.9) 



(6.10) 



with BY^^{^i) defined in (|7.1| ). Fb is the 5-meson decay constant and rjB is the QCD factor 
defined in (|6.8|). 

In the extensions of the SM with minimal fiavour violation (MFV) and without contributions 
from new operators it is useful to generalize ( |6.9| ) to 



/i^2 _ r2 



487r2 



M^mBF'sBB{V:,VuYTlB^Fu 



where 



'IB 



Vb 



Qncw 



(6.11) 



(6.12) 



Here r]]^^ and Sq™ describe new physics contributions in analogy to rjB and So{xt), respectively. 
That is 

«(5)(/i6) 



-6/23 



1 + 



An 



-J, 



(6.13) 



where C^^ifJ't) is the new physics contribution to CY^^ifJ't) and ASq"^^ results from new physics 
contributions without QCD corrections. 

In more complicated models in which new fiavour-violating couplings are present and the 
full set of operators ( |2.1[ ) is relevant, it appears to be most useful to evaluate new physics 
contributions using simply 



G'p 
16^ 



(6.14) 



with Ci{fi) evaluated by means of the formulae in sections and |^. Similar comments apply to 
the — system with obvious replacements. 
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7 Phenomenological Applications 

7.1 Hadronic Matrix Elements for — Mixing 

The matrix elements {K^\Qi{^)\K^) = {Qi{fi)) contributing to — mixing can be written 

as 



/2(/i)mxF|5fL(/i) 



24 
1 



where 



R{^^) 



rriK 



(7.1) 
(7.2) 
(7.3) 
(7.4) 
(7.5) 

(7.6) 



^rrisin) + main) J 

and Fx is the J'T- meson decay constant. Let us calculate the non-perturbative parameters 
B^{ji) using the lattice results of discussed in 0. In the Landau RI scheme (LRI) the 
Bf{fj,) factors are given by 



5[^(/.) 

'B!^\^^) 



LRI 
LRI 
LRI 



[^i(/^)]lri' 

[^5(^)]lRI> 
[^2(/i)]LRI, 



Bf\^.) 



LRI 



LRI 



[^4(/i)]LRI> 

^-BM - '^-BM 



(7.7) 



LRI 



where Bi{fi), i = 1, . . . ,5 are the non-perturbative factors entering the matrix elements in the 
operator basis of [^. 

In order to find the matrix elements ( [7.1|) -( [73D in the NDR scheme we use the results of , 
which allow us to relate the Bi factors in the LRI and NDR schemes. We find 

a(^)(/i). 



Bl^{^^) 



NDR 



NDR 



B^^i^,) 



LRI 47r 
LRI 477 



„VLL 



i?V-(/^) 



LRI 



LRI 



(7.8) 
(7.9) 
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NDR 



NDR 



NDR 



Bf\^.) 



LRI 



+ 



LRI 



+ 



LRI 



+ 



in 
An 



LRI 



12 



LRI 



LRI 



where 



,VLL 



-LR 



sSLL 



A^Lm-.NDR = 0.8785, 



= Ar 



-1.1288 -6.7726 \ 
LRi^NDR - Q_3069 10.8712 J ' 

5.6438 0.2140 \ 



SLL 

LRI^NDR 



.12.9387 2.6892;' 
Now, the Bi factors presented in [1^] read for /i = = 2 GeV as follows: 

[^i]lri = 0.60 ±0.06, 
[^3]lri = 1.05 ±0.12, [B 



[fi2]LRi = 0.66 ±0.04, 



4JlRI 



1.03 ±0.06, 



\B 



sJlri 



0.73 ±0.10. 



Using the central values for these parameters, we find by means of ( [7.7| 



^VLL 
^LR 



SLL 



LRI 
LRI 
LRI 



0.60, 
0.73, 

0.66, 



B^^ 



LRI 



B: 



SLL 



LRI 



1.03, 



0.40. 



(7.10) 
(7.11) 
(7.12) 

(7.13) 
(7.14) 

(7.15) 



(7.16) 



(7.17) 



Finally, setting a(^)(2 GeV) = 0.306 and using ( |7.8| )-( [7?T^ ) we obtain in the NDR scheme for 
^ = ^^ = 2 GeV 



^VLL 
^LR 
^SLL 



NDR 
NDR 
NDR 



0.61, 
0.96, 

0.76, 



NDR 



B. 



SLL 



NDR 



1.30, 



0.51. 



(7.18) 



We observe that the scheme dependence in the LR and SLL sectors is large, amounting to 
a shift of the Bi factors by roughly 30% and 20%, respectively. The corresponding scheme 
dependence in the VLL sector amounts to 2%. 
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Setting Fk = 160 MeV and tuk = 498 MeV we obtain at /i = 2 GeV 



(Qi^)ndr 



(Q2^)nDR 



(Qi'"'")ndr 



(Q2^'")nDR 



m,(2GeV) +md(2GeV)_ 

1 2 



0.26 ■ 10-^GeV^ 
-3.83 

7.77 

-3.79 

-6.10 



llSMeV 



llSMeV 



m,(2GeV) +mrf(2GeV) 

llSMeV 
m,(2GeV) +md(2GeV) 

llSMeV 
m,(2GeV) + md(2GeV) 



• 10"2GeV^ 
■ 10~2GeV^ 

2 

■ lO-^GeV^, 

2 

■ lO-^GeVl 



(7.19) 
(7.20) 

(7.21) 

(7.22) 

(7.23) 



7.2 AMk, AMb and £k 

Next we would like to present general formulae for the mass differences AM^ and AMb in the 
— and — systems and for the CP- violating parameter e^- In the case of AM^ 
and Ek our formulae are valid for the contributions of heavy internal particles with masses 
higher than M\y- The known SM contributions from internal charm quark exchanges and 
mixed charm-top exchanges []T1| have to be added separately. 
We have 



AMk = 2Re{K^\H^J^=^\K°) 



AMb = 2|(5°|i7,^^=2|s0\ 



\/2AMx 

The matrix element {K°\H^^=^\K'^) can be written as follows 



(7.24) 
(7.25) 
(7.26) 



487r2" 



M^ruKFl P,VLL(Cr^(/i,) + Cr"(/it)) 



VRR/ 



+PfL(Cf^(/Xi) 



CiS^^(/Xi)) + PfL(CfL(/x,) 



^CT\l^t))\ (7.27) 
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with 



VLL 



pLR 



jLR 



eff 



eff 



PI 



jSLL 



8 



eff 2 

3 

eff ~ 2 



eff 



eff 



(7.28) 

(7.29) 
(7.30) 
(7.31) 
(7.32) 



In the case of the SM and MFV models one can use ( |6.9| ) and ( |6.11| ) instead of ( [7.28| ). In 
writing these formulae we have absorbed the CKM factors into Ci{nt)- The effective parameters 
[Bi{fiL)]es are defined by 



[BKf^L)] 



rriK 



eff 



BK^^l 



18.75 



115 MeV 



n 2 



BKi^l). (7.33) 



In the case of — mixing one has to make the replacements f^b and rriKF^ msF^. 

Then in the case of B^ — B^ system 

"4.4 GeV 



mif^b)] 



TUB 



eff 



1.44 



1 2 



B^fib), (7.34) 



^rubifib) + mdifib) 
with an analogous formula for the B^ — B^ system. 

We would like to emphasize that these formulae together with the QCD factors rjij presented 
in Section 3 are valid for any extension of the SM. In particular the coefficients are universal. 
New physics contributions enter only the coefficients CKfit)- The latter have to be evaluated 
in the NDR scheme in order to cancel the scheme dependence of the universal coefficients . 
In the process of multiplying Cf^Ht) and terms 0{a1) have to be removed. 

It is instructive to evaluate the coefficients P" for the — system. Setting fii = 2 GeV, 
^Ws ^ "^^^ MeV, nisifiL) + fndifii) = 115 MeV and using the values for Bf in ( [7.18| ) we find 



(7.35) 
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pLR = _3g 1^ pLR ^ 3^ ^7_3g^) 

pSLL ^ _i3 1^ pSLL ^ _32 2. (7.37) 

We observe that the coefficients P^"^ and P^^^ are by two orders of magnitude larger than P^^^. 
This originates in the strong enhancement of the QCD factors rjij for the LR and SLL (SRR) 
operators and in the chiral enhancement of their matrix elements as seen in ( [7.33| ). Consequently 
even small new physics contributions to Cl"^{fit) and Cf^^{fit) can play an important role in 
the phenomenology H, |^ . 

In the case of -B" — B'^ mixing the chiral enhancement of the hadronic matrix elements in the 
LR and SLL sectors is absent. Moreover, the QCD factors rjij are smaller than in the — 
mixing. Consequently the coefficients P}'^ and pf^^ are smaller in this case. As lattice results 
are not yet available for all the relevant hadronic matrix elements in the B system [|I5| we will 
set Bi{mb) = I. Taking uib = 5.28 GeV, fXb = GeV, mfe(/ife) + ma{[ih) = 4:A GeV and 
= 226 MeV we find 

pVLL _ Q_g4^ (7 3g) 

pLR ^ -i_Q2, P^^ = 2.46, (7.39) 
pSLL ^ _i_47^ pSLL ^ _2 (7 40) 

8 Summary 

We have presented analytic formulae for the QCD renormalization group factors relating the 
Wilson coefficients Ci{fit) and Ci{fi), with fit = 0{mt) and fi < fit, of the AF = 2 dimension 
six four-quark operators Qi. The formulae presented in section 3 are given in the NDR scheme 
but are otherwise universal and apply to the Standard Model and all its possible extensions. 
In order to complete the evaluation of AF = 2 amplitudes, the QCD factors presented here 
have to be combined with the Wilson coefficients Ci{fit) evaluated in a given model at the short 
distance scale fit and with the hadronic matrix elements {Qi{fj)) evaluated at ft = fif,, fi = fii 
or fi = fix dependently on the process considered. Ci{fit) and {Qi{fi)) have to be computed 
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in the NDR scheme in order to obtain a renormahzation scheme independent answer for the 
physical amphtudes. 

We have also presented analytic formulae for the QCD factors relating the matrix elements 
{Qi{2 GeV)) and {Qi^Hk)) with hk < 2 GeV. These formulae allow the comparison of the 
matrix elements obtained in lattice simulations with those obtained in approaches which use 
lower renormahzation scales. 

Our numerical analysis shows that the renormalization-group effects are very large in the 
LR and SLL sectors. This in particular is the case for the elements [?72i(/i)]LR, ['722(Ai)]LR, 
[?7ii(/i)]sLL, ['7i2(/u)]sLL and [?722(a*)]sll and is in accordance with the previous analyses [|, P]. 
The NLO corrections amount typically to 5-15% except for the elements [ri2i{lj)]hK, ['722(At)]LR 
and [7722(Ai)]sLL, where in the case of /i = IGeV they can reach 25%. As a result of this pattern 
the renormahzation group effects strongly enhance the Wilson coefficients and Cf^^ and 
strongly suppress Cf^^ with respect to their values at Ht- The corresponding effects in C^^^ 
and are substantially smaller. 

Finally we have presented expressions for the mass differences AMk and AMb and the 
CP-violating parameter ex in terms of the non-perturbative parameters Bf and the Wilson 
coefficients Ci{^t)- These formulae include renormahzation group effects at the NLO level and 
allow to calculate straightforwardly AMk, AMb and ex in any extension of the SM once the 
Wilson coefficients Ci{fit) and the non-perturbative parameters Bf are known in the NDR 
scheme. In the case of — mixing we have presented the results for [i?f (2 GeV)]p^j-)p^ using 
the lattice results obtained in the LRI scheme. The corresponding results for the B system 
should be available this year. 
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Appendix 



A One-Loop and Two-Loop Anomalous Dimension Ma- 
trices 



We give below the one-loop and two-loop anomalous dimension matrices. The two-loop expres- 
sions are given in the NDR scheme 0. 

^(l)VLL ^ _7 + 



.(O)VLL 



-fy-'— = 4, 



1 



(0)LR 



2 




12 
-16 



7 



(1)LR 



71 _ 22 J 
3 9 J 

^-2/ 



198 - f/ 

1343 _|_ 68 



(A.l^ 



7 



(O)SLL 



-10 

-40 



34 I ' 
3 / 



1459 I 74 



7 



(l)SLL 



_35 ]_ f 

36 54 



6332 I 584 f 
" 9 9 



2065 394 J 
9 27 



B Expansion of the evolution matrices U in a. 



Recall the renormalization group equation to which (|2.12| ) is the formal solution. At NLO it 
reads, written in terms of a^, 



2(3q a. 



+ 



2/5, 



27^' -—7^ 



1 

2A) 



1 

47r 



f/(/il,/i2). 



(B.2) 



At the leading order, where only the term oc l/a^ is kept, ( p.2|) has the (exact) solution 



At the next-to-leading order one has |l| 



exp 



7 



(0)T 



2/3o 



log 



(B.3) 



(B.4) 



where higher orders in the parentheses have been omitted. The algorithm for constructing the 
matrix Jj can be found in |l|, H, |13|. Eq. ( p.4|) holds in a theory with a given number / of 
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active quark flavours. When evolving across a quark threshold, as for instance in ( p.l4| ), one 
finds 



f/(/iL,/it) 



1 + 



47r 



47r 



1 + 



An 



(0) 
/=5 



In light of the fact that 0{a1) and higher terms have been dropped in ( |B.2|) and 
adopt the convention 

( as(/i2)\^ 



log(r^) = 



(B.5) 

), we 

(B.6) 
(B.7) 



and drop all 0{a1) and higher terms in ( p.5|) and similar expressions. For the desired two- and 
three-step evolution matrices, one obtains 



(0) 



r(0) 



r(0) r7(o) 



(B.8) 



(0) 



47r 



isf/fi^}! + V3U% {I - J3) Ui% - VsV^uflsUi^l] , (B.9) 



(0) 



r(0) 



K 



Aw 



:b.io) 



where we have suppressed some obvious arguments in the LO evolution matrices U'^^^ in order 
not to unnecessarily clutter the expressions. 

C The Evolution Matrix ?7(/it, /is) 

For completeness we give here the elements of the evolution matrix U {fit, f^s) in a / = 6 flavour 
theory with fig > fit- The renormalization group evolution from fig down to fit can even be 
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included as in ( |1.4| ) when fis is only by a factor of two higher than mj. However, it is only 
necessary when fis > 4:mt in order to avoid large logarithms. 

The formulae given below are not as general as the ones given in section ^ They apply 
only to the evolution of new physics contributions which do not involve SM particles except 
for the number of quark flavours entering ag and the anomalous dimensions of the operators 
This is for instance the case considered in p, |^ in which squarks and gluinos have been 
integrated out at a scale fig ^ yUt- On the other hand the renormalization group analysis of 
charged Higgs contributions with Mh± ^ irit would be more complicated as both and top 
can be simultaneously exchanged in box diagrams. Integrating out first and subsequently 
the top would introduce bilocal structures for fit < fi < fig quite analogous to the study of 
charm contributions to — 1^ mixing [0, |l^ . We find then 
VLL-Sector 



VLL 



VLL 



LR-Sector 



Vnif^t) 



SLL-Sector 



LR 
LR 
LR 
LR 
LR 
LR 
LR 
LR 



6/21 

1.3707(1 - r/6)r/, 



6/21 



3/21 



-24/21\ 



24/21 



(C.ll) 
(C.12) 



(C.13) 
(C.14) 
(C.15) 
(C.16) 

0.9219 rjQ-^^/^^ + r]Q^'^^ (-2.2194 + 1.2975 r/e) , (C.17) 
1.3828 (7/6^^/21 (C.18) 
(-10.1463 + 0.8650 7/6) +r/6"^^/^^ (-6.4603 + 15.7415 r/e) , (C.19) 
0.9219 ?76^^/2^ + r?6~^^^^^ (9.6904 - 10.6122 r/g) • (C.20) 



^11 Vt 



SLL 



1.0153r76-o•69l6-0.0153r/°■^«6^ 



(C.2i: 
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SLL 




(C.22) 




SLL 


V/D 10 / ' 


(C.23) 


Tl'nn { UA 
'122 xr'tj 


SLL 


= 1.0153r;°-'^6^ - 0.01537/6-°-'''', 


(C.24) 




SLL 


= riQ-^-^^^^ (5.6478 - 6.0350 T^g) + T^g"'^^^^ (0.3272 + 0.0600 T^g) , 


(C.25) 




SLL 


= 7?6"°-"" (10.7494- 37.9209776) + r/6°-^^" (41.2556 - 14.0841 r^g) 


,(C.26) 




SLL 


= 776°-^^" (-0.0618 - 0.0315 7^6) + r/6-°-"" (0.0454 + 0.0479 7^6) , 


(C.27) 




SLL 


= 776"°-^^^^ (0.0865 + 0.3007 r/g) + r/g^'^^^^ (-7.7870 + 7.3999 r/g) . 


(C.28) 



Here fit = 0{mt) and 77g = af^fXs) / Oi^^KfJ't) ■ These results together with those presented in 
section 3 and 4 allow to find f/(/x, /is) with < fit, see (p..4|). 
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